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Integrated quantum photonic waveguide circuits are a promising approach to realizing future pho- 
tonic quantum technologies. Here, we present an integrated photonic quantum technology platform 
utilising the silicon-on-insulator material system, where quantum interference and the manipula- 
tion of quantum states of light are demonstrated in components orders of magnitude smaller than 
previous implementations. Two-photon quantum interference is presented in a multi-mode interfer- 
ence coupler, and manipulation of entanglement is demonstrated in a Mach-Zehnder interferometer, 
opening the way to an all-silicon photonic quantum technology platform. 



Introduction 

Quantum information technologies offer completely 
new approaches to encoding, processing and transmit- 
ting information. By harnessing the properties of quan- 
tum mechanics, such as superposition and entanglement, 
it has been shown possible to realise fundamentally new 
modes of computation [UH], simulation [3, 4| and com- 
munication ,5|, as well as enhanced measurements and 
sensing 6 . Of the many prospective physical systems 
in which to encode quantum information, photons are 
a particularly promising approach due to their proper- 
ties of low noise, easy manipulation and low transmission 
losses. To date, quantum photonics integrated circuits 
have been realised in low index contrast waveguide ma- 
terial systems, such as silica 15] and silicon-oxy-nitride 
[9] • Such technologies offer benefits in terms of low prop- 
agation losses, but their associated large bend radii limits 
the scalability and usefulness of this technology. 

Here we present silicon quantum photonic waveguide 
circuits utilising the silicon-on-insulator material sys- 
tem, where quantum interference and the manipulation 
of quantum states of light were demonstrated in com- 
ponents orders of magnitude smaller than previous im- 
plementations. Quantum interference of indistinguish- 
able photons was realised in multi-mode interference cou- 
plers, and manipulation of multiphoton entanglement 
was demonstrated in a Mach-Zehnder interferometer. 



Integrated quantum photonics 

Traditionally, photonic quantum information experi- 
ments have been implemented using bulk optical ele- 



ments, with photons propagating in free space. Although 
many proof-of-principle experiments have been reported, 
such an approach rapidly becomes impractical as the 
complexity of the quantum optical circuits increases, 
thereby making them inherently unscalable and confin- 
ing them to research laboratory optical tables. In addi- 
tion, sub- wavelength stability is critical for reliable oper- 
ation of many quantum circuits due to the necessity of 
often complex networks of nested interferometers. How- 
ever, developments over the past few years have overcome 
these bottlenecks through implementation of integrated 
quantum circuits, allowing quantum information science 
experiments to be realised that are inherently stable and 
orders of magnitude smaller than their equivalent bulk 
optic implementations. Integrated quantum photonic cir- 
cuits have been demonstrated in a host of different ma- 
terial systems, including silica-on-silicon [H] [TO], direct 
write silica (7J [IT] , silicon oxy-nitride [9] , lithium niobate 
[12] and gallium nitride [13], all realized in low index- 
contrast waveguide structures. Two-photon quantum in- 
terference has been reported in both silica and gallium 
nitride [101 113] , whilst manipulation of quantum states of 
light has been demonstrated using silica-on-silicon [14] 
and lithium niobate |12j ; the later being also a promising 
material for photon pair generation |15j . The silicon-on- 
insulator wire waveguide technologies offer a further level 
of control, stability and miniaturisation, and also allow 
routes for on-chip generation and detection of photons. 
The high refractive index contrast provides compactness 
unrivalled by any other photonic material system, allow- 
ing a dramatic reduction in the footprint of quantum 
circuits and the integration of complex circuits on a sin- 
gle chip. Thermo-optic heaters, p-n junction modulators 
and induced x 2 nonlinearities |16] allow for dynamically 
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reconfigurable circuits with the possibility of fast phase 
control. The high third-order nonlinearity of silicon en- 
ables on-chip generation of quantum states of light via 
spontaneous four wave mixing |17j . thus allowing inte- 
gration of single photon sources and waveguide circuits 
on the same chip. Recently, high efficiency on-chip sin- 
gle photon detectors integrated with silicon waveguides 
have been reported [18] . Demonstration of quantum in- 
terference and manipulation of quantum states of light 
in a silicon waveguide circuit is the next critical step to 
realizing photonic quantum technologies in silicon. In this 
work we report quantum interference and phase manip- 
ulation of one and two photon states using thermo-optic 
phase shifters and multimode interference devices on a 
silicon chip, opening up the way to a CMOS compati- 
ble silicon-based photonic quantum technology platform, 
where sources, detectors and circuits can all be realized 
monolithically on the same integrated chip for applica- 
tions in communication, computing and metrology. 



a) b) 




FIG. 1: MMI and MZI a) SEM image of a multimode 
interference device, b) Schematic diagram of waveguide circuit 
with voltage controlled phase shifter, c) Illustration of the 
cross-section of the single-mode waveguide. 



Experimental detail 

The silicon integrated waveguide circuits were fabri- 
cated from a silicon-on-insulator (SOI) substrate with 
a silicon thickness of 220nm. The single-mode optical 
waveguides had a width of 450 nm, and a top cladding 
of silicon dioxide. Input and output coupling to lcnscd 
fiber was achieved using spot-size converters (SSC) com- 
prising a 300 /im-long inverse-taper with a 200 nm tip 
width and a 4x4 [im 2 polymer waveguide. The silicon 
waveguiding structures were defined by 248 nm lithogra- 
phy and formed by dry etch processing. Typical propa- 
gation losses of 3.5 dB/cm and SSC losses of 2 dB/facet 
where observed, leading to device losses in the range of 
6 — 10 dB. To realise the beam-splitter-like operation re- 
quired for quantum interference, a 2x2 multi-mode in- 
terference coupler was designed (using FDTD simula- 
tions) and implement with dimensions of 2.8 /^m x 27 
fj,m and input tapers of length 3 nm and width 1 /im 
(Fig. Qa)). Mach-Zchnder interferometers were formed 
from two multimode interference couplers and a 200 /im- 
long thermo-optic phase shifter, shown in Fig. |TJb) . By 
varying the voltage across the phase shifter it was possi- 
ble to tune the internal phase delay of the Mach-Zehnder 
interferometer. Pairs of indistinguishable single photons 
at a wavelength of 1550 nm were generated from a type-I 
spontaneous parametric down conversion source (SPDC). 
A bismuth borate BiB 3 Og (BiBO) nonlinear crystal was 
pumped by a 775 nm, 150 fs pulsed Ti-Sapphire laser, 
and degenerate photon pairs were collected from two di- 
ametrically opposite points on the SPDC cone using sin- 
gle mode fibers, and coupled into the chip using lensed 
fibers. One arm of the collection system was on a motor- 
ized delay stage to provide a tunable delay between the 
photon pair. The 1550nm photons were coupled out of 



the chip using single-mode lensed fiber and detected us- 
ing two superconducting single photon detectors (SSPD), 
having system detector efficiencies of 5% and 15% |19ll2T)] . 
Coincidence detection between the two output ports of 
the device was performed using a custom made (FPGA 
based) counting logic with a 5 ns coincidence window. 



Quantum interference in a 2x2 multimode 
interference coupler 

Quantum interference lies at the heart of linear quan- 
tum photonics and is fundamental to the implementa- 
tion of any photonic quantum technology, describing at 
it's most basic level the interaction of indistinguishable 
photons incident at a beamsplitter. In the case of two 
identical photons, each entering separate ports of a 50:50 
beamsplitter device, the amplitude probabilities of both 
photons exiting different ports are zero — giving a su- 
perposition state of either both photons exiting one out- 
put port, or both photons exiting the other output port. 
This quantum interference occurs for any physical imple- 
mentation of a two-port beamsplitter device, and for in- 
tegrated waveguide circuit two typical implementations 
are the directional coupler and multimode interference 
coupler. 

Multimode interference (MMI) devices are based on 
the self-imaging principle, and compared with directional 
couplers they have improved tolerances to fabrication, a 
wider spectral bandwidth of operation, and can be ex- 
tended to realise NxN multiport devices with several in- 
put and output ports [3T] . To demonstrate quantum op- 
eration of the silicon-on-insulator 2x2 MMI coupler, in- 
distinguishable photon pairs were launched into each of 
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the two input waveguides. At the two output ports, the 
two-photon coincidence counts were monitored. By vary- 
ing the temporal delay between the two input photons, 
the indistinguishability of the input state could be con- 
trolled and the Hong-Ou-Mandel (HOM) dip [22 could 
be observed (see inset of Fig. [2]) . A maximum dip visibil- 
ity of 80% ± 3% was observed, demonstrating quantum 
interference. 

A limited visibility of 80% was achieved in part due 
to multi-photon events in the source, in part due to in- 
trinsic losses within the MMI device, and also due to 
residual distinguishability of the generated photon pairs. 
The indistinguishability of the SPDC generated photon 
pairs was characterised by performing a HOM experi- 
ment with a bulk-optic beamsplitter, giving a maximum 
visibility of 95.5% ± 1.5% for low pump power. Due to 
the chip insertion losses (a total of 9 dB including in- 
put/output coupling) and detector efficiencies (5% and 
15%), it was necessary to pump the SPDC source in a 
regime where high count rates could be achieved, which 
gave rise to multi-photon events that reduced the visi- 
bility of the dip. In order to quantify the effect of these 
multiphoton events, the visibility at different pump pow- 
ers was measured (see Fig.j2|, clearly demonstrating im- 
proved visibility for reduced photon-pairs per pulse. By 
fitting to a model which accounts for multi-pair gener- 
ation and system losses (see Appendix), an extrapolate 
nominal visibility of 88% ± 3% at low pump power was 
demonstrated. 

The 7.5% discrepancy between the measured beam- 
splitter visibility of 95.5% ± 1.5% and the nominal MMI 
visibility of 88% ±3% can be accounted for by the intrin- 
sic losses within the MMI device. These losses do not en- 
force a ip = 7r phase shift between the probability ampli- 
tudes for both photons being reflected and both photons 
being transmitted. A perfectly balanced 2x2 MMI should 
operate in exactly the same way as a perfect beamsplitter 
with the following scattering matrix: 



1 


' i i 


i 


' i i 


71 


i -i 


~V5 


1 



In the presence of losses, we account for the loss modes 
by embedding the scattering matrix of the MMI in a 
larger 4x4 unitary matrix written as: 



J 

(rja) 1 
(rja) 1 
(rja) 1 



(rja) 2 
e iv (rja) 1 - 
e w (rjaf- 

■ R 1 

e 1 ' (rja) 1 



(rja) 2 
ie (rjay- 



(rja) 2 

" R - 

e 11 (rja) "• 



where < 77 < 1 is the reflectivity of the MMI, < 
a < 1 is the loss parameter, rj = 1 — 77, a = 1 — a, ip is 
the internal phase of the MMI, and j3 are the phases 
of the loss modes. Consequently, the phase between the 
probability amplitude for both photons being reflected 
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FIG. 2: Visibility as a function of the probability of 
getting a photon pair Visibility of the Hong-Ou-Mandel 
experiment is plotted as a function of the probability per pulse 
of generating a photon pair. Solid line is a theoretical fit 
using a model which accounts for multi-photon events and 
losses before and after the MMI. Inset shows the two-photon 
interference plot for the highest visibility measurement. 



and both photons being transmitted is not necessary it 
anymore. The unitary condition gives a boundary on the 
internal phase ip as a function of the loss a: 
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The observed nominal visibility of 88% could be ex- 
plained by an intrinsic MMI loss greater than 0.8 dB, 
which could yield ip = 2.74 instead of tt phase shift. In- 
deed, FDTD simulations predict an intrinsic device loss 
of 0.5 dB, which is within experimental error of the 0.8 dB 
calculated by the model. By reducing the intrinsic loss 
of the MMI to 0.2 dB (through improved device design 
and fabrication) , the lower bound on the visibility would 
improve to 99.5%. It is worth mentioning that mode 
propagation simulation shows that the distorsion of the 
single photon wavepacket (~100 /im coherence length in 
free space) is insignificant in this experiment and does 
not contribute to a reduction of the visibility. 



Two-photon entangled state manipulation 

The ability to actively prepare and measure arbitrary 
quantum states is crucial for the implementation of many 
quantum information processing experiments. Manipu- 
lation of a dual rail-encoded qubits — a single photon 
in an arbitrary superposition of two optical channels — 
requires control of the relative phase and amplitude be- 
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tween the two optical paths. Arbitrary single qubit op- 
erations can be implemented using just two Hadamard 
gates (beamsplitter) and phase shifters, with the sim- 
ple Mach-Zehnder Interferometer (MZI) considered to be 
most basic fundamental building block required to re- 
alise any arbitary N-mode linear quantum photonic cir- 
cuit [25] , 

To demonstrate the operation of this fundamental 
building block, a silicon-on-insulator MZI was formed 
from two 2x2 MMI couplers and a thermo-optic phase 
shifter (see Fig. [ljb)). The MZI on it's own can apply 
only a subset of single qubit operations, but any arbitrary 
single qubit operation can be achieved by simply adding 
another phase-shifter element before and after the MZI 
[Hj. This device was characterised in both the single- 
photon and two-photon regime. 

In the first instance, single photons where input into 
port a (see Fig. Wb) ), and as the internal phase of 
the MZI was changed, by applying a voltage across the 
thermo-optic phase shifter, the probability of detect- 
ing photons at the output port c varied sinusoidally as 
P c = i[l — cos(cf))], with a periodicity of 2tt. The ob- 
servation of this classical interference fringe represents 
the ability to transform a single-photon input in mode a, 
into a superposition state across modes c and d, given by 
the transformation 1 10) -> cos(4>/2) 1 10) + sin{<j>/2) |01), 
where (j> represents the phase shift within the interferom- 
eter. In the second instance, the two-photon |11) state 
was used, with one photon injected into port a and one 
into port b. This input state leads to quantum inter- 
ference at the first MMI couplers and, after propagating 
through the phase shifter, transforms to the two-photon 
entangled state ^(|2) |0) + exp 2 ^ |0) |2)). Quantum in- 
terference at the second MMI coupler enables analysis 
of this quantum state, giving the probability output of 
coincidence detections P c -d — |[1 — cos (20)], with a pe- 
riodicity of 7r — half that of the single-photon case. 

The measured single-photon and two-photon fringes 
are shown in Fig. [3] for an applied voltage change from 
V to 4.5 V, resolving 1.5 fringes for the single-photon 
case, and 3 fringes for the two-photon case. The two- 
photon fringe had a visibility of V = 81.8%± 1.3%, which 
is greater than the threshold V t h = 1/V2 required to 
beat the standard quantum limit, demonstrating quan- 
tum metrology and the ability to achieve sub shot noise 
limited measurement. 



Discussion 

This is the first demonstration of interference and ma- 
nipulation of quantum states of light in silicon integrated 
quantum circuits, and is a fundamental step for further 
miniaturisation of photonic quantum circuits. All previ- 
ous integrated waveguide quantum circuits have relied on 
weakly guided waveguide structures, with a typical index 
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FIG. 3: Single photon and two-photon fringes a) Sin- 
gle photon count rate as a function of voltage applied to the 
thermo-optic phase-shifter for the input state |10). b) Two- 
photon coincidence count rate as a function of applied volt- 
age for the input state |11). The x axis is the applied voltage 
across the thermo-optic phase-shifter. The top figure show 
the single photon fringe. The y axis represents the number 
of single photon recorded at one output of the MZI. The bot- 
tom figure shows the two-photon fringe. The y axis represents 
the number of coincidence counts obtained for the two-photon 
fringe. 



contrast for the silica-on-silicon waveguide technologies of 
0.5%, resulting in a bend radius of greater than 10 mm 
[14] . This low index contrast and large bend radius re- 
sult in physically large circuits and components, making 
future implementations of complex quantum circuit ar- 
chitectures impractical. Due to the high index contrast 
of the silicon-on-insulator material systems, it is possible 
to achieve bend radii of below lOfim - three orders of 
magnitude smaller than the silica material system. Indi- 
vidual components can be dramatically reduced in size, 
such as the 2x2 MMI coupler presented here which is 
40 times smaller in length (27/xm) compared to 1.1 mm 
length of the equivalent device in silica [23] . This enables 
circuits of significantly greater complexity to be realised 
whilst still maintaining a small chip size, but potentially 
at the cost of higher waveguide losses. 

The key functions required to realise a fully integrated 
linear-optic quantum technology are the generation, de- 
tection, interference and manipulation of quantum states 
of light, all on a single optical on-chip. Generation and 
detection of single photon states in silicon wire waveg- 
uide devices has been previously demonstrated [TTJ, [TB] , 
and this work demonstrates the feasibility of on-chip in- 
terference and manipulation of quantum states of light in 
silicon wire waveguide circuits. The high x 3 nonlinearity 
of silicon and the high modal confinement of the silicon 
wire waveguide enables efficient photon pair generation 
via spontaneous four wave mixing, with a recent demon- 
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stration showing indistinguishable photon pair genera- 
tion from two independent silicon wire waveguides with 
an external HOM dip visibility of 73% [25]. High effi- 
ciency detection of single-photons has also recently been 
demonstrated in silicon wire waveguide incorporating su- 
perconducting nanowire detectors, with an internal de- 
tector efficiency of 94% [TH]. These recent demonstra- 
tions along with the work presented here represent the 
basic building blocks required to realise an integrated 
photonic quantum technology platform where quantum 
states can be generated, manipulated, interfered and de- 
tected all on the same circuits, opening up new possibili- 
ties in quantum information science and applications. Ul- 
timately this CMOS compatible technology could be in- 
tegrated with conventional microelectronic circuits, pro- 
viding on-chip driver circuits and fast logic. 



Conclusion 

Quantum interference and manipulation of quantum 
states of light in silicon wire waveguide circuits has 
been demonstrated. A maximum HOM dip visibility 
of 80% ± 3% was observed for a 2x2 multimode in- 
tereference coupler, and a nominal visibility in the ab- 
sence of multi-photon terms of 88% ± 3% was calcu- 
lated. Internal MMI losses were shown to be the be 
dominant mechanism for the reduced nominal visibil- 
ity. An integrated Mach-Zehnder interferometer with 
a thermal phase shifter demonstrated on-chip entangled 
state manipulation with a two photon fringe visibility of 
81.8% ± 1.3%. These results pave the way for the reali- 
sation of fully integrated photonic quantum technologies 
in silicon. 
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APPENDIX 

This appendix describes the various details of the 
model used to give the fit in Fig. [2] 

Multiphoton treatment of quantum interference 

The state produced by the SPDC source is assumed to 
be a non degenerate squeezed state and can be written 
as: 

oo 

|*) = vT^]Tf»|n> 

n=Q 

with < £ < 1 where £ is the squeezing parameter, 
n is the number of photon pairs generated, |n) |n) is the 
number of photons in the two spatial modes of the source. 
We label Pci,Pc2 and Pec the probabilities of detect- 
ing per pulse, respectively a single photon at detector 1, 
a single photon at detector 2, and a coincidental event 
between the two detectors CI and C2. 



Those probabilities can be expressed as: 



Pa = # 7 aEe 2n (i-(i-')in 

oo 

p C 2 = v/FaE^f 1 -' 1 -*)") 

oo 

Pec = v / FaE^( 1 '( 1 ^i)")( 1 -( 1 ^)") 

with r\\ and r\i being the overall efficiencies (includ- 
ing collection, transmission and detection) of respectively 
channels 1 and 2. 

We measured the quantities CI (number of detection 
events per second at the detector 1), C2 (number of detec- 
tion events per second at the detector 2) and CC (number 
of coincidental detection events per second at the detec- 
tor 1 and 2) for different pump powers. 

From those datasets, we ran a minimisation methods to 
relate the different count rates CI, C2, CC obtained from 
the pump powers Ik to the channel efficiencies rji and r/ 2 , 
and the (, 2 {Ik) parameter using the following equations: 

C1(I) = fxPci(I) 
C2(I) = fx P C2 (I) 
CC(I) = f x P cc (I) 

where / = 80 MHz is the repetition rate of the pulsed 
laser and / is the intensity of the pump. 

We then launched the state |>F) for different pump pow- 
ers into the MMI coupler, and recorded quantum interfer- 
ence patterns with visibilities V(I). In order to infer the 
visibility at low intensity, we applied the model described 
below. 



Theoretical visibility 

The theoretical visibility is expressed as 1 — -p^ where 
Pj is the probability to get a coincidental event when 
the photons generated in the two modes are indistin- 
guishable, and Pjj is the probability to get a coincidental 
event when the two photons are distinguishable. If we 
used a true single photon source in an ideal lossless cir- 
cuit, Pi = 0, Pd — 0.5 and V = 1. However, we used a 
squeezed state |\E') as the input state, and the setup ex- 
hibits collection losses, propagation losses and non uni- 
tary detection efficiencies. We modelled those losses by 
adding a virtual beamsplitter at each input and output 
port of the MMI coupler with reflectivity's r\A,f]B and 
Vg'iVd as shown in Fig. [4j And we trace over the loss 
modes La, Lb, Lc, Pd- AH that is needed is to compute 
Pi and Pd taking into account these losses. 
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FIG. 4: Input and output losses modeling of a 2x2 MMI cou- 
pler 

We start with a pure squeezed state with a density 
matrix p = |^} ($|, we first propagate the state through 
the two input beamsplitters, we then trace over the loss 
modes L^and Lb- We further propagate the obtained 
reduced density matrix through the MMI coupler, as- 
suming a 50/50 splitting ratio. And we finally propagate 
the output state through the last two loss beam splitter 
and trace over the modes Lq and Ljj. 

By writing pi ou t as the final density matrix, we com- 
pute the probability to get a coincidence by calculating: 

oo oo 

P. i = J~] ^ («| (m\ piout \n) \m) 

n—1 m—1 

We implemented an algorithm to compute this proba- 
bility to an arbitrary precision and found that computing 
the terms up to n, m — 9 was accurate for the different 
pump power we used. Adding higher order terms did not 
changed significantly the results. 



For the case where the photons are distinguishable, 
we start with \^dist) defined below and assume that the 
photons are orthogonal. In practice, one arm is delayed 
with respect to the other and we consider the photons 
launched in mode a at time t\ and the photons launched 
in mode b at time ta, giving: 

oo 

Then applying the same evolution as previously for the 
indistinguishable photon case, we compute the probabil- 
ity Pjj of getting a coincidence when the photons are 
distinguishable . 



Fit of the experimental data 

We recorded the visibility as a function of the pump 
power. For each intensity, we compute the associated 
squeezing parameter (£ ). The collection and detector 
efficiencies and the device losses are known. The only 
free parameter over which we need to minimize is the 
visibility. We therefore use an equation of the form: 

, )= aPj(0 + (l-a)P D (0 
^' Pd (£) 

where a is the parameter which quantifies the over- 
lap between the photons which includes residual spectral 
entanglement from the source and imperfections in the 
MMI coupler. Using this model, we extrapolate the nom- 
inal visibility that would be expected at very low pump 
power. 



